Molecular CP-violating magnetic moment 
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A concept of CP-violating (T,P-odd) permanent molecular magnetic moments fi CP is introduced. 
We relate the moments to the electric dipole moment of electron (eEDM) and estimate /i CP for 
several diamagnetic polar molecules. The moments exhibit a steep, Z 5 , scaling with the nuclear 
charge Z of the heavier molecular constituent. A measurement of the CP-violating magnetization 
of a polarized sample of heavy molecules may improve the present limit on eEDM by several orders 
of magnitude. 
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It is a common knowledge, that heteronuclear diatomic 
molecules possess a static electric dipole moment aligned 
with the internuclear axis n, D = Dn. For a dia- 
magnetic molecule, however, there is no similar mag- 
netic moment. As demonstrated below, an existence 
of such a magnetic moment would violate both parity- 
transformation (P) and time-reversal (T) discrete sym- 
metries. Because of the compelling CPT theorem, an 
observation of this magnetic moment would provide an 
evidence for the CP violation Q, Q ■ CP-violation, al- 
though observed in particle physics, still remains a mys- 
tery, as much stronger CP-violating mechanisms may be 
required to explain the matter-antimatter asymmetry of 
the Universe. 

Here we introduce the molecular CP-violating mag- 
netic moments, /x CP = ^ CP n. We propose a mea- 
surement of n via detection of ultra-weak magnetic 
fields generated by a polarized sample of diamagnetic 
molecules. For several molecules we evaluate /i CP and 
express them in terms of the permanent electric dipole 
moment (EDM) of electron, d e . A measurement of non- 
vanishing molecular CP-violating moments would reveal 
the elusive electron EDM (eEDM) (EDMs violate both 
T and P symmetries). While no EDMs have been found 
so far, it is anticipated that the next generation of ex- 
periments may finally discover the EDMs. Indeed, most 
supersymmetric extensions of the Standard Model of ele- 
mentary particles predict eEDMs that are within a reach 
of planned and on-going experimental searches (see, e.g., 
a popular review 3]). The present limit on d e comes from 
an atomic Tl beam experiment H, 



d e (Tl) < 1.6 x 10 



-27, 



(1) 



Here we propose an experimental search for the CP- 
violating magnetic moments of heavy polar molecules. 
We argue that the limit on fi CP derived from such ex- 
periments would imply constraints on d e that are several 
orders of magnitude better than the present limit • I n 
principle, the experiments can be carried out with any 
diamagnetic polar molecules. However, there is no par- 
ticular advantage in using polyatomics, and we restrict 
our consideration to polar diatomic molecules. 



In the remainder of this paper, unless specified oth- 
erwise, we use atomic units |e| = h = me = 1 and the 
Gaussian system for electro-magnetic equations. In these 
units, the Bohr magneton is fiB = a/2, where a « 1/137 
is the fine structure constant, and the unit of magnetic 
field is me 2 e 5 /?i 4 w 1.72 x 10 7 Gauss. 

General considerations. Diatomic molecule is charac- 
terized by the projection £1 = (J-n) of the total electronic 
angular momentum J = L + S on the internuclear axis 
h. For a molecular state with a definite fl, the molecular 
magnetic moment is directed along h and, phenomeno- 
logically, we may construct the following combinations of 
the two vectors 



^ CP n 



HbG\\{J ■ n)h, 



(2) 



where /x CP and G\\ are numbers. For the Hund's case 



(a) G-factor is given by an expression G||f2 w A + 2S, 
where A = (L ■ h) and £ = (S ■ n) 5]. While the sec- 
ond term in is T,P-even, the /i CP n term violates both 
time-reversal (and thus CP) and parity. Indeed, under 
the time reversal the magnetic moment acquires a mi- 
nus sign, while h is T-invariant. Similarly, under parity 
transformation, fi is not affected, while h flips direction. 

Given a complete set of molecular states \k) (with en- 
ergies Ek), the magnetic moment fi CF of a state |0) can 
be computed as 



CP 



(M ■n)\k)(k\V^\0) 

En — Eh 



(3) 



where M is the operator of magnetic dipole moment, and 
the CP-violation is due to a T,P-odd interaction V CP . 
Here we consider eEDM as a source of CP-violation, so 
that V CP = -d e {jo - 1)70752 • £int, where 7* are the 
conventional Dirac matrices |l|. The electric field £j nt 
exerted upon the EDM is the strongest at the nucleus, 
leading to £i„t ~ (Z/r 2 ) r, where Z is the nuclear charge 
and r is the radius- vector of the electron with respect to 
the nucleus. The matrix element of V CP scales as Z 3 , so 
that the interaction V can be considered as localized at 
the heavier nucleus. Below we will evaluate the molecular 
sum © using an approach similar to the LCAO method 
(linear combination of atomic orbitals). 
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Note that Eq. © is expressed in the body-frame of the 
molecule. After ^ CP is found one has to average Eq. (J2J 
over rotations. In the external electric field (n) ^ and 
we get magnetization in the direction of the electric field. 
The second T,P-even term in Eq. J5J) does not contribute 
to this magnetization. For diamagnetic molecules (J • 
n) = Q = and this term vanishes. For paramagnetic 
molecules in the absence of the magnetic the levels with 
different signs of f2 are equally populated and this term 
is averaged to zero. 

We would like to detect a magnetization of a sample 
of polarized molecules due to CP-violating magnetic mo- 
ments. In this regard, it is beneficial to work with dia- 
magnetic molecules 57 — 0, so that the traditional T,P- 
conserving magnetic moment (the last term in Eq. @) 
does not contribute to the magnetization. Some of the 
molecules may still have non-zero nuclear magnetic mo- 
ments. However, the magnetization due to the nuclear 
moments in a macroscopic sample will average out to 
zero, since the nuclear spins have equal probabilities 
of orienting parallel or anti-parallel to the direction of 
internuclear axis. Another advantage of diamagnetic 
molecules is that they are also chemically stable allowing 
for higher sample densities and thus for a larger sample 
magnetization. 

To illustrate our qualitative approach to evaluat- 
ing CP- violating magnetic moments, consider a polar 
molecule CsF in its ground ^ state. Halides exhibit a 
chemical bond of a strong ionic character, and we model 
the CsF molecule as the Cs + ion perturbed by the elec- 
tric field £ of negative ion F~. The perturbing field at 
the Cs + is £ ps q/R%, where R e is the internuclear sepa- 
ration and q = 1 is the valency of Cs. The CP-violation 
is enhanced near the heavier atom and we may evaluate 
the magnetic moment as 

^ p (CsF)«/3 CP (Cs+)-|, (4) 

where CP {Cs + ) is a so-called CP-violating polarizabil- 
ity [10 of the Cs+ ion. Thus the molecular two-center 
problem is reduced to computing a one-center prop- 
erty — CP-violating polarizability of the heavier con- 
stituent. If both constituents of the diatomic molecule 
AB have comparable nuclear charges, then /z CP (AB) w 
[/3 cp (A(+«))-/? CP (B(-«))] qR- 2 , where q is the observed 
valency of the atoms. 

CP-violating polarizabilities. The effect of CP- 
violation on electromagnetic properties of the media 
has been considered by a number of authors (see, e.g., 
BUSES)- In particular, the atomic CP-violating (T,P- 
odd) polarizability /3 CP relates induced atomic magnetic 
moment /z at to the externally applied electric field £, 



For a spherically-symmetric system, the CP-odd polar- 
izability is a scalar quantity, i.e., the induced magnetic 
moment is directed along the applied E-field. 



In Ref. |?J we have computed and explored the atomic 
CP-violating polarizabilities for rare-gas atoms in their 
respective Sq ground states. Since the calculation re- 
ported here builds upon that work, let us briefly reca- 
pitulate our approach and main results. We related (3 
to the eEDM through third-order perturbation theory: 
the CP-odd polarizability of an atomic state <E>o can be 
represented as 

/3 CP = -2(<5>o\V CP KM z KD z \<5> ) 

-2(<5> \M Z KV CP KD Z \<1> ) (6) 

-2(<5>o\M z KD z KV CP \<5> ), 

where the resolvent operator 1Z = (e — Ha.t)~ , £0 be- 
ing the energy of the state <&o and H at being the atomic 
Hamiltonian. We have evaluated Eq.© using the Dirac- 
Hartree-Fock (DHF) approximation. Numerical evalu- 
ation has been carried out using a B-spline basis set 
technique 111]]- 

In Ref. |3| , we demonstrated that (3 GP has a doubly rel- 
ativistic origin: relativistic effects are essential for non- 
vanishing matrix elements of V (Schiff theorem) and 
also due to peculiar properties of magnetic-dipole opera- 
tor M. This doubly relativistic nature leads to a steep, 
Z 5 R(Z), scaling of (3 CP with the nuclear charge, where 
slowly-varying R(Z) is a relativistic enhancement factor 
(see Ref. Q for details). This result, together with the 
estimate (0J , immediately provides us with the Z-scaling 
of the CP-violating magnetic moments of diamagnetic 
molecules 

M CP cx Z 5 R{Z) , 

where Z is the nuclear charge of the heavier constituent. 
This Z-scaling is more substantial than the usual Z 3 R(Z) 
scaling ^3] of the effects of eEDM on the energy levels 
of paramagnetic molecules (radicals). As in the tradi- 
tional EDM searches, it is beneficial to search for non- 
vanishing CP-odd magnetic moments with molecules in- 
volving heavy atoms. Diamagnetic contribution JSJ to 
CP-odd polarizability is universal. However, for the sys- 
tems with unpaired electrons this contribution is masked 
by a larger paramagnetic term which is linked to the 
electronic angular momentum @ . The advantage of the 
diamagnetic systems is the much lower magnetic noise, 
which may be crucial for the experiment. On the theo- 
retical side, calculation of fj, for diamagnetic molecules 
may be simpler and more reliable than for para- or fer- 
romagnetic crystals. 

Results of calculations. In Table [i] we present the re- 
sults of our calculations of CP-violating magnetic mo- 
ments for several diatomics: CsF, BaO, TIF, PbO, and 
BiF. These diamagnetic molecules possess the ground 
state. The heavier atoms of these diatomic pairs are met- 
als, and we assume that the molecules exhibit a pure case 
of ionic bond, i.e., these heavier atoms fully lend their 
valence electrons to their electronegative companions (F 
and O) and become closed-shelled ^SO ions. The second 
and third columns of the Table list the resulting heavy 
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atomic ions with their nuclear charges, and in the fourth 
column we present our computed values of CP-violating 
polarizabilities of these ions. Finally, we combine ionic 
(3 CP with the equilibrium internuclear separations (see 
Eq.@) and obtain an estimate for the molecular CP- 
violating magnetic moments. Our sign convention in ex- 
pression /x CP = fi CP h is such that the unit vector h is 
directed from the heavier to the lighter nucleus. Notice 
that we express the /i CP in terms of eEDM. Similarly, 
/i CP can be expressed in terms of the T,P-odd electron- 
nucleon couplings; a simple prescription of how to relate 
the present results to the strength of such couplings is 
given in Ref. 0. 



Molecule Ion 


ZIon /3 UP 


(Ion)/d e 


H /de 




CsF Cs+ 


55 


-3.0[-2] 


1.5[- 


•3] 


BaO Ba ++ 


56 


-2.3[-2] 


3.4[- 


■3] 


TIF T1+ 


81 


2.9[-l] 


-1.9[- 


-2] 


PbO Pb++ 


82 


3.2[-l] 


-4.9[- 


-2] 


BiF Bi+ 


83 


4.8 


-3.2[- 


-1] 



TABLE I: Molecular CP-violating magnetic moments, 
u CP /d e , divided by the eEDM for several diamagnetic 
molecules. The values of /i CP /d e are dimensionless, d e and 
/^ CP being expressed in the Gaussian atomic units. The sec- 
ond, third, and the fourth columns list the heavier ion in the 
molecule, its nuclear charge and its CP-violating polarizabil- 
ity, /3 CP /d e , in units of fi 4 /(me 2 e 5 ). Notation x[y] stands for 
x x 10 y . 



For all the considered molecules, the internuclear sepa- 
ration R e w 2A, and thus the internal molecular fields are 
comparable. More significant is the effect of increasing 
CP-violating polarizabilities (the fourth column of Ta- 
ble P) as one progresses to heavier elements. This trend 
is largely due to the Z 5 scaling of f3 . Yet, there is an 
order of magnitude of difference between /3 CP for Pb ++ 
(Z = 82) and Bi + (Z = 83). A part of this large enhance- 
ment lies in a softer excitation spectrum of Bi + and thus 
smaller energy denominators in Eq. JBJ. Also, while solv- 
ing the DHF equations we assumed that the outer shell 
of Bi + ion has the 6p 2 / 2 electronic configuration. In gen- 
eral, however, the ground state of Bi + would contain a 
combination of Qpy 2 an d 6P3/2 configurations. Since the 
P1/2 states couple to EDM strongly, while p 3 / 2 orbitals 
contribute at a much smaller level, we expect that our 
result for (3 CP of Bi + is somewhat overestimated. 

Results of Table [I] should be considered as a qualitative 
estimate for another reason as well. The expressions (0J 
and @) are based on atomic wavefunctions | <!>,), instead 
of the molecular wavefunctions of the defining expres- 
sion ©• An underlying assumption is that the molecular 
wavefunctions \i) in the vicinity of the heavier atomic ion 
can be expressed perturbatively as 



where |\&o) is the wavefunction of the lighter ion (we 
left out excitations from |^o) as being non-essential for 
computing /i CP ). Certainly, ab initio molecular-structure 
calculations of CP-violating magnetic moments are desir- 
able to depart from this simple perturbative picture. To 
motivate more sophisticated calculations, let us evalu- 
ate a feasibility of an experimental determination of ji 
based on our qualitative estimates. 

Proposed experiments. The molecular CP-odd mag- 
netic moments are tiny. For example, combining the 
present limit on the eEDM (JIJ with the computed value 
of jjF p /d e , we obtain for BiF, 

^ CP (BiF) < 2.4 x lCT 37 crg/Gauss. (7) 

While this is a remarkably small value, only 2.6 x 10~ 17 
of the electron magnetic moment, measuring such small 
magnetic moments seems to be possible with the modern 
magnetometry. 

The CP-violating magnetic moments can be deter- 
mined by measuring the ultraweak magnetic field gen- 
erated by a sample of molecules. Because of rotations, 
the body-fixed /x CP moment averages to zero in the lab- 
oratory frame. Experimentally, one needs to apply a po- 
larizing electric field £ po \ to orient the molecules along 
the field. For the efficient polarization, the coupling to 
the field must be stronger than the rotational spacing, 
D£ P oi > 2.B, where B is the rotational constant. For the 
ground vibrational state of BiF, the rotational constant is 
B sa 0.231 cm" 1 , requiring the application of the polariz- 
ing E-field of a few kV/cm. Another consideration comes 
from thermal averaging over rotational levels [l3| , and it 
is beneficial to work at very low temperatures to have the 
smallest possible number of populated rotational levels. 

Suppose we have a spherical cell with polarized dia- 
magnetic molecules; the permanent CP-violating mag- 
netic moments of the molecules produce a macroscopic 
magnetization of the sample generating an ultraweak 
magnetic field. The maximum value of the magnetic field 
at the surface of the cell can be expressed as 



where n is the number density of the sample. One could 
measure this induced magnetic field and set the limits on 
the eEDM or other CP-violating mechanisms. Clearly, 
one should increase the number density to enhance the 
signal. However, condensing polar molecules with ionic 
bonds leads to a crystallization of the sample. To main- 
tain the individuality of the molecules, one could em- 
ploy low-temperature matrices of rare- gas atoms with 
molecules embedded inside the matrix |l4j . The ma- 
trix isolation is a well established technique in chemical 
physics. For chemically stable molecules, the number of 
guest molecules per host atom (matrix ratio) , could be as 
high as 1/10 0], i- e -j with the techniques of matrix isola- 
tion, one could attain the number densities of molecules 
in the order of 10 21 cm -3 . 



4 



With the sample number density of 10 21 cm -3 and 
/i CP of BiF derived from the present limit on eEDM, 
from Eq.JHJ we obtain a generated B-field of B ~ 2 x 
10~ 15 Gauss. This ultraweak field can be measured 
within a month of integration time at the present best 
sensitivity limit p{ of 5 X 1CT 12 G/\/Hz. A relatively 
small 0.3 cm 3 -volume sample has been used in that ex- 
periment. The projected theoretical limit |16( of this 
method is 10~ 13 G/a/Hz. More optimistic projected sen- 
sitivity of 3 x 1CT 15 G/VHz is given in Ref. [l7|. With 
this projected sensitivity, we find that the present limit 
on the eEDM Q may be improved by three orders of 
magnitude, 

d e (BiF, projected) «2x 1(T 30 e • cm. (9) 

In Ref. Q, we have considered a CP-odd magnetiza- 
tion of a sample of liquid Xe caused by an externally 
applied electric field J^J. Here we focused on a similar 
magnetization due to molecular CP-odd magnetic mo- 
ments and we find that the BiF experiment has a sub- 
stantially better sensitivity to eEDM. This enhancement 
is due to (i) larger nuclear charge of Bi (Z = 83) than 
that of Xe (Z = 54) and (ii) much larger E-field applied 
to heavy atom/ion: in case of BiF, the internal molecu- 
lar field is ~ 4 x 10 8 V/cm, while in liquid Xe the E-field 
is limited by the breakdown strength of 4 x 10 5 V/cm. 
This large difference in the maximum attainable labo- 



ratory field and the internal molecular field 18] is ex- 
ploited in more conventional searches for EDMs with 
molecules It is worth emphasizing that 

the experiment considered here is based on a bulk mag- 
netization of diamagnetic molecules, while the conven- 
tional searches for eEDM determine energy splittings in 
individual paramagnetic molecules. 

To summarize, we introduced a concept of molecular 
CP-violating (T,P-odd) magnetic moments, /i CP . We 
related these magnetic moments to eEDM and we esti- 
mated fi CF for a number of diamagnetic polar molecules. 
We demonstrated that /j, cp exhibit a strong Z 5 depen- 
dence on the nuclear charge of the heavier molecular con- 
stituent. Finally, we evaluated a feasibility of setting a 
limit on the eEDM by measuring ultra-weak magnetic 
fields produced by a polarized sample of diamagnetic 
molecules. We found that such an experiment may im- 
prove the present limit on the eEDM by several orders of 
magnitude. 
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